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Abstract 

We present a calculation of the NLO QCD corrections to heavy flavor photo- 
production with longitudinally polarized beams. We apply our results to study 
the spin asymmetry for total charm quark production which will be used for a 
first direct determination of Ag by the COMPASS experiment. We also briefly 
discuss the main theoretical uncertainties inherent in this calculation. 
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1 Introduction 

Despite significant progress in the field of spin-dependent EjlS-the polarized 
gluon density Ag remains almost completely unconstrained uau by presently 
available DIS data. An important role plays here the lack of any direct con- 
straint on Ag from other processes. Upcoming spin experiments will thus put 
a special emphasis on exclusive measurements to provide further invaluable 
information for a more restrictive analysis of polarized parton densities in the 
future. In this context heavy quark (Q — c, b) photoproduction is considered to 
be one of the best options to pin down Ag because in LO only the photon-gluon 
fusion (PGF) processu (an arrow denotes a longitudinally polarized particle) 

19 -> QQ (1) 

contributes. Since LO estimates of (Q) are rather unreliable, and since we 
already know from the unpolarized NLO calculation ErQ that the corrections 
are sizeable and that the 'clean picture' of (Q) is obscured by new, light quark 
induced NLO subprocesses, the knowledge of the polarized NLO corrections 
is mandatory for a meaningful extraction of Ag. In what follows wc will briefly 
highlight the main steps and results of our calculation 0. 



2 Technical Framework 

The NLO QCD corrections to the PGF mechanism in consist of three 
parts: (i) the one-loop virtual corrections, (ii) the real corrections 75 — > QQg, 
and (iii) a new genuine NLO production mechanism jq (q) — > QQq (q). In the 
calculation of (i)-(iii) one encounters UV, IR, and mass (M) singularities which 
are removedJpy renormalization, in the sum of (i) and (ii), and by factorization, 
respectively 13. To make all these singularities manifest we choose to work in 
the framework of n-dimensional regularization. 
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The required polarized squared matrix elements for ([!]) and (i)-(iii) 



AlMl 2 = - 
1 1 2 



\M\ 2 (++) - \M\ 2 (+- 



(2) 



where the ± denote the helicities of the incoming particles, are obtained by 
projecting onto the helicity states of the bosons (photons_pr giuons) and quarks 
using the e M „ p(T tensor and the 75 matrix, respectively 13. By taking the sum 
insteac.£jf the difference in (|^) we fully agree with the known urmolarized 
resultsErO and the abelian, 77 — > QQ, part of our results agrees withB as well. 

The presence of 75 and t^ vpa in the polarized calculation introduces some 
complications since these objects have no unique continuation to n 7^ 4. In the 
HVBM prescription, which we use, the usual n-dim. scalar products k ■ p are 
accompanied by their [n— 4)-dim. subspace counterparts k ■ p ('hat momenta'). 
These terms deserve special attention when performing the 2 — > 3 phase space 
integrations. For single-inclusive heavy quark production, which we consider 
here, only a single hat momenta combination p 2 shows up in A \Mf and the 
appropriately modified 2 — > 3 phase space formula schematically readsQ 

1 rl x (n-6)/2 

dPS 3 = dPS 3 ,„„p(0i, 9 2 ) x / dx (3) 

where x = 4(s4 + m 2 )p 2 / (s 2 sin 2 61 sin 2 62), B is the Beta function, m denotes 
the heavy quark mass, and S4 = s + ti -\-u±. 9ia are introduced to parametrize 
the momenta of the two not^bserved partonsQ. However, due to the appear- 
ance of m in x it turns out that all contributions due to p 2 are at least of 
0{n — 4) and hence drop out when the limit n — > 4 is takenQ. Tha remaining 
phase space integrationd then proceeds as in the unpolarized cased. 

Finally, it should be recalled that in the factorization procedure for (iii) 
one has to introduce the parton content of the polarized photon a which is 
experimentally completely unknown so far. A scheme independent result in 
0(a 2 a) can thus only be obtained for the sum of the 'direct' and 'resolved' 
photon contributions. The NLO corrections fox-the latter are unknown and, 
for the time being, have to be estimated in LOEj. 



"This differs from a calculation involving only massless particles. 
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Figure 1: The LO and NLO (MS) gluonic scaling functions as a function of r) = s/im? — 1. 

3 Numerical Results and Phenomenological Aspects 

The total photon-parton cross section in NLO can be expressed in terms of 
scaling functions [i = <?, q, q) 



Aa i7 (s,TO 2 ,/i/) 



aa s 

<m 2 



A/1%) 



A/1%) 



where A/^ 0) and A/f ] 



(4) 

A/j 7 stand for the LO and NLO corrections, respec- 
tively, fif denotes the factorization scale (for simplicity we choose (j, r — ///), 
and rj = s/4m 2 — 1. The scaling functions can be further decomposed depending 
on the electric charge of the heavy and light quarks, eQ and e q , respectively: 



A/ 97 (r?) 



e Q Ac g7 (?7) 



A/, 7 (t?) 



e Q Ac 97 (ry) 



3- Ad 



,(„) (5) 



with corresponding expressions for the A/ i7 . 

In Fig. 1 we present Acg 7 \ Acg 7 , and Acg 7 as a function of r\ in the MS 
scheme. For the discussion below it is important to notice that Acgty (solid 
line) changes sign at rj ~ 3. Upon adding the NLO terms, multiplied by a factor 
Aira s (see @)), the zero is shifted towards rj ~ 1. We also note that for 77 < 0.1 
the 0(a s ) corrections dominate over the LO result when we include that factor. 
Comparing our polarized results with the unpolarized ones (see Fig. 5 inB), 
one observes that for rj — > Ac ff7 — ► £37? 

implying that |M ff7 | (H — ) — ► 0. On 
the contrary, for r\ — > 00 the unpolarized NLO coefficients approach a plateau 
value □ dominating over the LO result while all polarized coefficients tend to 
zero here, implying that |M S7 | 2 (++) — > |M_ 7 | 2 (H — ). The numerically less 
important quark coefficients can be found in 
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Figure 2: j4^ p in LO and NLO for three sets of polarized parton densities <r^ p in (6) 
was calculated using the GRV 11 densities. The dotted line shows A^ p in LO using a NLO 
Ag (see text). The bar shows the uncertainty for such a measurement at COMPASS 12 . 

Using (|J) it is easy to calculate the experimentally relevant spin asymmetry 
for the total hadronic heavy flavor photoproduction cross section 

A^ p (S Jp ,m 2 ,Hf) = A<jQ p (S lp ,m 2 , n f )/<jQ p (S lp ,m 2 , fi f ) (6) 

as a function of the photon-proton c.m.s. energy S lp and where 




l 

Aa^ p (S 7P ,m 2 ,fi f ) = J dxAa fl (xS 7P ,m 2 ,(i f )Af p (x,(i 2 f ) (7) 

(with a similar expression for the unpolarized cross section a^ p ). 

In Fig. 2 we show the charm asymmetry A c lp in LO and. NLO, using m = 
1.5 GeV, fif = 2m, and three_sets of polarized densities Era, in the energy 
region relevant for COMPASS B (they will operate at w 10 GeV). The NLO 
corrections are large, depend strongly on S lp and do not cancel in the ratio 
(^) as one jnay naively expect. However, their origin is readily explained 13: 
For GRSV0 and GS (A) I and ^/K^ > 12 GeV the corrections stem from the 
shift of the zero in the gluonic coefficient function. A c lp changes sign at some 
J S 1P > 20 GeV and large NLO corrections in the vicinity of a zero are natural. 
For y/ 'S 1P < 12 GeV, where one probes Ag at x > 0.1, the corrections are, on 
the other hand, entirely due to the badly constrained Ag, more precisely, due 
to the large differences between the LO and NLO Ag in the two sets Era. Shis 
is illustrated by the dotted curve where we use the NLO GRSV gluon tl to 
calculate the LO A c lp . In this energy region the observed large corrections 

4 



20 



t> 



- lOGeV -__ 

I I I I | I I I I | I I I I | I I I I | I I I I | I I I I | I I I I | I I I I 



VS m = 100 GeV 

TP 



NLO 

- - LO 



\4 [GeV 2 ] 

Figure 3: The factorization scale dependence for the LO and NLO A<7^ p as defined in (7) 
using the LO and NLO GRSV densities 1 , respectively. 

thus should not be taken too literally. For the GS (C) seti the situation is 
more involved since here Ag oscillates as well. 

Finally we briefly discuss the importance of the main theoretical uncer- 
tainties. Fig. 3 shows the dependence of Aa^ p on the choice of fif. The 
improved scale dependence in NLO clearly underlines the importance of our 
NLO calculation 0. Light quark induced subprocesses contribute about 5% 
at ^S 1P ~ 10 GeV, but more for GS (C), and should be subtracted before 

Stracting Ag. The 'resolved' photon contribution was shown to be small in 
. More important is the unknown value of m, leading to shifts of about 30% 
when m is varied by 0.2 GeV around 1.5 GeV used in our calculations. Finally, 
for a reliable extraction of Ag from A^ p our knowledge of the unpolarized gluon 
has to be-improved as well since at large x the uncertainty in g(x, fj, 2 ) is rather 
sizeable 113. 
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